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Abstract 

Lovelock gravity is an important extension of General Relativity that pro- 
vides a promising framework to study curvature corrections to the Ein- 
stein action, while avoiding ghosts and keeping second order field equa- 
tions. This paper derives the greybody factors for D-dimensional black 
holes arising in a theory with a Gauss-Bonnet curvature-squared term. 
These factors describe the non-trivial coupling between black holes and 
quantum fields during the evaporation process: they can be used both 
from a theoretical viewpoint to investigate the intricate spacetime struc- 
ture around such a black hole, and for phenomenological purposes in the 
framework of braneworld models with a low Planck scale. We derive exact 
spectra for the emission of scalar, fermion and gauge fields emitted on the 
brane, and for scalar fields emitted in the bulk, and demonstrate how the 
Gauss-Bonnet term can change the bulk-to-brane emission rates ratio in 
favour of the bulk channel in particular frequency regimes. 



PACS Numbers: 04.70.Dy, 04.50.+h 



1 Introduction 



In any attempt to perturbatively quantize gravity as a field theory, higher- 
derivative interactions must be included in the action. Such terms also arise in 
the effective low-energy action of string theories. Furthermore, it has been argued 
that higher- derivative gravity theories are intrinsically attractive as in many cases 
they display features of renormalizability and asymptotic freedom. Among such 
approaches, Lovelock gravity [T] is especially interesting as the resulting equa- 
tions of motion contain no more than second derivatives of the metric, and as 
they are free of ghosts when expanding around flat space. Among the dimension- 
ally extended Euler densities used in the Lovelock Lagrangian, it is reasonable to 
consider that the four-derivative Gauss-Bonnet term is the dominant correction 
to the Einstein-Hilbert action [2] . The action therefore reads as : 

S GB = J d D xV=9 [R + a(R^ pa R^ - m^BT + R 2 )} , 

where a is a coupling constant of dimension (length) 2 , and G the D— dimensional 
Newton's constant defined as G = 1/M®~ 2 in terms of the fundamental Planck 
scale M*. In four dimensions, the Gauss-Bonnet term is a total derivative, but 
it becomes dynamical for D > 4. The thermodynamical properties of such black 
holes have been already studied both in asymptotically flat [3] and curved [3] 
spacetimes. Unfortunately, the temperature and entropy which have been de- 
rived are not sufficient to describe the detailed evaporation spectrum through the 
emission of Hawking radiation : the exact coupling between the black hole and the 
quantum fields, i.e. the greybody factors, must be accurately computed. To date, 
the greybody factors have been obtained for D-dimensional Schwarzschild 00 Ej 
(for a detailed review, see [S]), Reissner- Nordstrom jU], Kerr (TU1 EU El UH] and 
Schwarzschild-de-Sitter jHj black holes, for various types of emitted fields. 

Furthermore, it has been recently pointed out that small black holes could be 
formed at next-generation colliders during trans-planckian particle collisions ^3] 
if the Planck scale is of order a TeV, as is the case in some scenarios postulat- 
ing the existence of extra dimensions JB]. This idea has driven a considerable 
amount of interest in the framework of collider physics |T7], as well as in cosmic- 
ray physics ^H] as the same phenomenon could also occur during ultrahigh-energy 
neutrino interactions in the atmosphere of the earth. Most works consider those 
black holes as described by higher- dimensional generalizations of 4-dimensional 
line-elements that follow in the context of Einstein's linear theory of gravity. 
Studying the experimental consequences of the Gauss-Bonnet term, when in- 
cluded in the D-dimensional action, is a very challenging path strongly connected 
to experimentally probing quantum gravity. With arguments based solely on the 
dynamics of the temperature, it has been pointed out that for realistic values of 
the Gauss-Bonnet coupling constant and a low enough Planck scale, the LHC 
should be able to discriminate between a pure Einstein action and a second or- 
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der Lovelock action upon the successful detection of the Hawking radiation 
emission spectrum from such a black hole. 

To confirm this stimulating hypothesis it is nevertheless necessary to calculate 
the exact Gauss-Bonnet greybody factors values that have not been computed 
to date. The aforementioned work ^Hj used an approximate expression for the 
greybody factors to compute the emission spectrum from a Schwarzschild-Gauss- 
Bonnet black hole; this approximate expression ignored altogether its dependence 
on the energy and spin of the emitted particle as well as its dependence on the 
Gauss-Bonnet coupling constant itself. In this work, we calculate through nu- 
merical analysis the exact values of the Gauss-Bonnet greybody factors and the 
corresponding emission spectra for scalar, fermion and gauge boson fields re- 
stricted to live on our 4-dimensional brane. The spectrum for the emission of 
scalar fields in the bulk will also be studied, and the bulk-to-brane ratio will 
be accurately calculated. The dependence of the spectra on the Gauss-Bonnet 
coupling constant will be studied in detail, and footprints of its value would be 
looked for in both bulk and brane emission channels. As we will see, the imple- 
mentation of the exact Gauss-Bonnet greybody factors will lead to modifications 
in the Hawking radiation emission spectra compared to the ones derived in the 
pure Schwarzschild case, and can also change the so-far derived results in the 
literature for the bulk-to-brane emission rates ratio. 

In what follows, we will work in the context of the scenario with Large Extra 
Dimensions [16] . and assume that the higher-dimensional spacetime is empty and 
thus flat. Also, the horizon value r# of the produced black holes will be assumed 
to be much smaller than the (common) size 1Z of the extra dimensions; under 
this assumption, the extra spacelike dimensions can be considered to be non- 
compact as the usual 3-dimensional ones. The energy of the emitted particle u> 
will be taken to be always much smaller than the black hole mass Mbh so that 
no significant back-reaction to the black hole background will take place after the 
emission of a particle. The frequency of the emitted particle must also be bounded 
from below, i.e. uj ^> l/7£; this guarantees that the wavelength of the particle is 
much smaller than the size of the extra dimensions and the particle can therefore 
be considered as a genuinely higher- dimensional one. Although the upper limit 
on the frequency of the black hole will be always ensured in our analysis, for the 
sake of presenting complete spectra, the lower bound will be ignored and restored 
at the end of our analysis. 

Section 2 of this paper sets the general framework for our analysis. The ab- 
sorption cross-sections, i.e. the greybody factors, are then numerically computed 
for particles on the brane and in the bulk in sections 3 and 4, respectively. Sec- 
tion 5 presents the final Hawking radiation emission spectra for both brane and 
bulk channels, and the bulk-to-brane emission rates ratios is explicitly given. We 
discuss our results and summarise our conclusions in section 6. The numerical 
values of the greybody factors produced in this analysis, that might be useful for 
experimental searches, can be downloaded from the URL given in [2Uj . 
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2 General Framework 



A higher-dimensional, neutral, static, spherically-symmetric black hole within the 
framework of a curvature-squared, Gauss-Bonnet gravity is described by the fol- 
lowing line element, first derived in an d then extended to a different topology 
in g]: 

ds 2 = h(r) dt 2 -^- r 2 dn 2 D _ 2 , (1) 
where dQ 2 D _ 2 is the line-element over a unitary (D — 2)-dimensional sphere and 



h(r) - k i r2 1±i l i 64^-3X0-4)6^ 

is the metric function. In the above, corresponds to the total (D — 2)- 

dimensional solid angle. In [21], it was shown that a Schwarzschild-like space- 
time within Gauss-Bonnet gravity exhibits two branches corresponding to the ± 
signs entering the metric function. The (+)-sign branch describes a Schwarzschild 
black hole with a negative mass parameter embedded in an Anti-de Sitter Universe 
whereas the (— )-sign branch corresponds to an asymptotically flat Schwarzschild 
black hole with a positive mass parameter, nevertheless both branches are char- 
acterized by a positive ADM energy, as it was shown in [22] • The authors of Ref. 
jl] extended the solution of to different topologies for the black hole's event 
horizon : for k = 1, the event horizon is elliptic whereas it becomes flat for k = 
and hyperbolic for k = —1 j2H|. In this work, we restrict our analysis to the 
(— )-sign branch with an elliptic topology for the event horizon, as it corresponds 
to the usual Schwarzschild black hole solution. The (+) branch will be hereafter 
ignored and the k parameter will be set equal to +1. 

As in higher-dimensional general relativity, Schwarzschild-Gauss-Bonnet black 
holes present one event horizon located at r# and defined as the solution of the 
algebraic equation h(r) = 0. The thermodynamical properties of those black 
holes have been extensively studied in jl]. For non- rotating, neutral black holes, 
the thermodynamics is described by three quantities, all parametrized by the 
horizon radius: the mass Mbh 

MBH ~ 8nGT(^l) l 1 + ^ J' (3) 

where T stands for the Euler Gamma function, the temperature Tbh given by 
the surface gravity evaluated at the event horizon radius, 

rp = (0-3) jr 2 H + a(D-5)(D-4)} 
BH 4irr H [r 2 H + 2a(D-4)(D-3)} ' 1 ' 



and the entropy Sbh, defined by 



i dM B H 

zbh = I • (5) 

J-BH 
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The Hawking radiation emission process by the black hole given in Eq. JT} 
takes place in the higher-dimensional space-time with a spectrum given by the 
generalization of the four- dimensional one obtained in 



dt ~ Zj< exp (u/T B h) - (-l) 2s (2tt) d - 1 ' 

In the above equation, j stands for the total angular momentum quantum number 
and s for the spin of the particle. In this work, we will focus on massless particles 
for which \k\ = u>; in this case, the expression for the flux is greatly simplified 
as the phase space term reduces to an integral over the energy of the emitted 
particle 

dtdu - ^e* V {u/T BH )-{-l) 2s (2n) D - 1 ' 

The function aj a (u), usually called "greybody factor" or alternatively "absorp- 
tion/emission cross section", when accurately taken into account, makes the re- 
sulting radiation spectrum more complex than the usual blackbody law. This 
additional factor arises because of the non-vanishing probability for the emitted 
particles to be reflected on the gravitational potential barrier and re-absorbed by 
the black hole. The cross section is computed within a full quantum mechani- 
cal framework to determine first the tunnelling transmission probability |-4j>| 2 
in the aforementioned Schwarzschild-Gauss-Bonnet background. The cross sec- 
tion for a given angular momentum quantum number is then deduced from the 
generalization of the optical theorem 



2°-^-^T {{D - 2)/2) (2j + D-3)(j+D- 4)! 2 
°-i» = X J\ ^ • (8) 

As it is clear from the above, the Hawking radiation emission spectrum is 
defined by two quantities, the black hole temperature and the particle's grey- 
body factor. From Eq. (J3J, we may see that the temperature has an explicit 
dependence on both the Gauss-Bonnet coupling constant and the dimensionality 
of spacetime. As we will see, the greybody factor depends, too, upon the underly- 
ing gravitational theory and thus encodes valuable information on both of these 
parameters; as expected, it also depends on properties of the emitted particle 
such as its energy and spin. This elaborate dependence of the greybody factor 
on a number of fundamental parameters of the theory is the main motivation 
for performing the exact calculation of its value and of the corresponding black 
hole radiation spectrum within the framework of a D-dimensional Gauss-Bonnet 
gravity. 
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3 Radiation on the brane 



3.1 General equation for particles on the brane 

Particles confined on the brane propagate into a background whose geometry is 
induced by the bulk curvature. This induced geometry on the 4— dimensional 
brane is simply given by fixing the values of the extra azimuthal angular coor- 
dinates, 9i = 7r/2, and leads to the projection of the D— dimensional metric on 
the 4— dimensional slice that describes our world. The projected metric has the 
form: 

dv^ 

ds 2 = h(r)dt 2 -^—-r 2 (d9 2 + sin 2 9dp 2 ) , (9) 
h{r) 

where the metric function h(r) remains the same as in the higher- dimensional line- 
element given in Eq. (JTJ). Consequently, the profile of the curvature along the 
three non-compact spatial dimensions contains a footprint of the bulk curvature. 

The equations of motion describing the behavior of particles with spin re- 
stricted to our 4— dimensional brane are deduced by using the Newman-Penrose 
formalism [2H1I2Z!- Under the assumption of minimal coupling to gravity, these 
equations can be combined into one master equation for particles with spin 
s = 0, 1/2, 1. In [H] (see also [TU] ) , this master equation was derived for particles 
confined on a brane embedded in a D— dimensional Kerr black hole background. 
In the limit of zero black hole angular momentum, the Kerr metric reduces to the 
Schwarzschild one, and the induced background on the brane assumes the form 
given in Eq. Using the following factorized ansatz for the propagating field: 

V 8 (t, r, 9, ip) = e~ iu)t e~ im{p ^~ s P s (r)Y^(9) (10) 

where A = hr 2 and j is the total angular momentum quantum number, the radial 
part of the master equation of motion is found to have the form: 

. n d / a1 ,dP s \ ( uj 2 r 2 n . isur 2 ti ~\ n 

In the above, A = j(j + 1) — s(s — 1) is derived from the master angular equation 
satisfied by the spin- weighted spherical harmonics Y™- 

The aforementioned master equation has been successfully used to derive the 
Hawking radiation emission spectrum on the brane for scalar, fermion and gauge 
boson fields by a D— dimensional Schwarzschild black hole 0013, and for scalar 
fields by a D— dimensional Schwarzschild-de Sitter black hole ^3]. Recently, its 
complete form, that includes a non-vanishing black hole angular momentum, was 
used to study the emission of scalar fields on the brane by a D— dimensional 
Kerr black hole fj^]- In our case, since the line-element of a D— dimensional 
Schwarzschild- Gauss-Bonnet black hole has a Schwarzschild-like form, Eq. 
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Figure 1: Gravitational potential seen by scalar particles propagating in the 
background of an induced Schwarzschild-Gauss-Bonnet black hole. Upper Panels: 
The number of dimensions is fixed at 5 and the Gauss-bonnet coupling constant 
at 0.1 M~ 2 on the left and at 10 M~ 2 on the right. Lower Panels: The number 
of dimensions is fixed at 10 and the Gauss-bonnet coupling constant at 0.1 M~ 2 
on the left and at 10 M~ 2 on the right. 



can be consistently used to study the emission of scalar 1 , fermion and gauge 
boson fields emitted by such a black hole on the brane. 

Writing the master equation (JTTj) for scalar particles in an alternative form 
allows to validate our approach to be followed in the next section. By using 
the tortoise coordinate r*, defined by h(r)dr* = dr, and redefining the unknown 
radial function by U(r) = rP (r), the scalar equation takes a Shrodinger-like 
form: 

d 2 U . . {£(£+ 1) h'( r y 



dr ^ ' Kr)[^ L + -^-)U = u/U, (12) 

where £ is the orbital angular momentum. From the above equation, the gravita- 
tional potential seen by scalar particles can be easily read, and it is depicted in 
Fig. ^for various values of the fundamental parameters of the theory. This poten- 

1 The equation of motion for a scalar field propagating in a D— dimensional Schwarzschild- 
Gauss-Bonnet black hole background was recently studied in to obtain the evolution of the 
scalar field in time and its quasinormal modes. 
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tial, for Schwarzschild-Gauss-Bonnet black holes, tends to zero both at the black 
hole's event horizon and at spatial infinity, due to the vanishing of the metric 
function in the first case and the asymptotic flatness of space-time in the second. 
The corresponding near-horizon and far-field asymptotic states can be therefore 
defined in terms of plane waves, a result that allows the use of the optical theorem 
for the computation of the absorption cross-section. Figure [0 clearly shows that 
the height of the potential barrier is increasing for modes with increasingly larger 
angular momentum number, and a similar effect arises as the dimensionality of 
spacetime increases, too. On the other hand, the height of the barrier is consider- 
ably lowered as the value of the Gauss-Bonnet coupling increases, with the effect 
being more significant for low values of D. 

3.2 Numerical results for the absorption cross-section 

Equation (fTT]) has been attacked analytically in the literature, and results have 
been produced for the absorption cross-section (see for instance Refs. [5^ IHl Hlj) 
under the assumption of low-energy emitted particles. As the gravitational back- 
ground, however, becomes more complicated, our ability to solve this equation 
analytically is greatly reduced. Moreover, the need to produce exact results, free 
from any approximations and valid throughout the energy regime, dictates the 
use of numerical analysis. The intricate shape of the metric function (j2J) in this 
particular case makes Eq. (jll)) impossible to solve analytically. Therefore, in this 
work, the calculation of the absorption cross section is performed numerically, for 
particles with spin s = 0, 1/2 and 1. 

The numerical integration of Eq. (jllj) starts from the black hole horizon, 
where appropriate boundary conditions are applied, and extends to spatial infin- 
ity. In order to deal with finite numbers and avoid divergences at the horizon, the 
numerical integration starts at = rn+s instead of r#, where e is small enough 
to induce a negligible error. As mentioned before, the asymptotic solution near 
the horizon is of a plane-wave form. In order to simplify the numerical analysis, 
the P s function and its derivative are fixed, at r ini , to the following values : 



which ensure that the no-outgoing modes boundary condition is satisfied, and fix 

2 

the flux at the black hole event horizon to \P S \ =1. 

The solution of the master equation is then propagated from the event horizon 
until a high enough value of the radial coordinate r, to be considered as spatial 
infinity. At this asymptotic regime, the wave function must be described by 
spherical waves [2] : 




(13) 



dP s 



(14) 



dr 





(15) 
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where B\ is the amplitude of the ingoing modes and B2 the amplitude of the 
outgoing modes. The numerical results are then compared to the above asymp- 
totic solution, and a fit is performed in order to extract the amplitudes B\ and 
B 2 . As it is obvious from the above asymptotic solution, for scalar particles with 
s = 0, the solution of the master equation at infinity contains both ingoing and 
outgoing modes. For particles with non-vanishing spin though, the complete so- 
lution follows by solving a set of two equations : one for the upper component 
s = +\s\, which is mainly composed of ingoing modes, and one for the lower 
component s = —\s\, which is mainly composed of outgoing modes [[HI- The 
absorption probability can be found through the ratio of the incoming flux at 
infinity over the one at the horizon, therefore only the solution for the ingoing 
modes is necessary for our study, and thus only the upper component of Eq. 
is solved. The absorption cross-section is then easily computed in terms of the 
absorption probability, through Eq. (jHJ) evaluated for D = 4 to account for brane 
emission : 



To obtain the final absorption cross-section, all the angular momentum contribu- 
tions have to be summed. Fortunately, as it is shown in Fig. ^ the gravitational 
potential barrier increases with j, and the high angular momentum contributions 
are strongly suppressed. The sum is therefore performed until rank j = m such 
that the error due to higher momentum contributions becomes irrelevant. 

In the case of particles with non-vanishing spin, another difficulty arises due 
to the very small value of the outgoing modes in the upper component of the 
wave function. As these modes are extremely small at the event horizon, a tiny 
error in the numerical solution can easily be mixed with an outgoing solution 
and can consequently contaminate the ingoing one. As a result, the numerical 
integration is not stable as it strongly depends on the boundary conditions. This 
contamination does not arise for fermions, but it becomes significant for gauge 
boson fields, and the absorption cross-section cannot be robustly computed by 
integrating Eq. A solution to overcome this problem has been proposed in 

[7j and consists in solving the equation of motion for a new unknown radial wave 
function F(y), defined as P\ = yF(y)e~ tu;r , with y = rjru- In terms of these 
new variables, the wave equation for gauge boson fields becomes : 



with the following boundary conditions having replaced the ones in Eqs. (|13|)- 
flUD : F(l) = 1 and dF/dy{\) = + l)/2ur H - 

3.2.1 Absorption cross-section for scalars 

The absorption probability is defined as the ratio of the ingoing flux JF at the 
event horizon over the one at infinity. For scalar particles, the flux is given by 




(16) 




(17) 
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Figure 2: Absorption cross-section for scalar fields on the brane for Schwarzschild- 
Gauss-Bonnet black holes as a function of the energy of the particle ujth- Left 
Panel : the Gauss-Bonnet coupling constant a is fixed at 0.1 M~ 2 and the number 
of dimensions D at {5, 6, 10}. Right Panel : the Gauss-Bonnet coupling constant 
takes the values {0, 0.1, 1, 10} M~ 2 and the number of dimensions is fixed at 6. 

the radial part of the conserved current J M = 2hr 2 Xm(^ qC^^o) , which leads to 
T = 2hr 2 Tm{l-%jj-r 2 Pl ). The absorption probability is found to have the form 




where, in the middle part, we have used the fact that the ongoing flux at the 
black hole horizon has been fixed to 1 by the boundary condition (I13|) . The last 
part of the above equation offers an alternative but equivalent way of calculating 
| A, | 2 through the reflection probability given in terms of the amplitudes of the 
ingoing and outgoing modes at infinity. 

Figure El shows the numerical results we have obtained for the absorption 
cross-section in units of nr 2 H as a function of the dimensionless parameter loth- 
The two panels depict results for : a fixed value of the Gauss-Bonnet coupling 
constant and variable D, the first one, and a fixed value of D and variable a, the 
second. In both panels, when urn 0, the absorption cross-section tends to a 
non- vanishing value (47rr H ), which does not depend on the number of dimensions 
or on the value of the Gauss-Bonnet coupling constant. This behaviour follows 
the same pattern seen in other cases of asymptotically-flat spherically-symmetric 
black hole backgrounds ED] emitting scalar radiation: in the low-energy 
regime, the absorption cross-section reduces to the area of the black hole horizon 
- for brane emission, this is simply ^mr 2 H . 

On the contrary, the behaviour of the cross-section in all other energy regimes 
strongly depends on both the number of dimensions and on the Gauss-Bonnet 
coupling constant : as for the case of Schwarzschild and Schwarzschild-de Sitter 
black holes, an increase in the number of dimensions leads to a decrease of the ab- 
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Figure 3: Absorption cross-section for fermion fields on the brane for Schwarzs- 
child- Gauss-Bonnet black holes as a function of the energy of the particle urn- 
Left Panel : the Gauss-Bonnet coupling constant a is fixed at 0.1 M~ 2 and the 
number of dimensions D at {5,6, 10}. Right Panel : the Gauss-Bonnet coupling 
constant takes the values {0, 0.1, 1, 10} M~ 2 and D is fixed at 6. 



sorption cross-section On the other hand, an increase of the Gauss-Bonnet 

coupling constant leads to an increase of the absorption cross- sect ion. Both of 
these results could be inferred from the behavior of the gravitational potential, 
that was derived in section 3.1, with respect to D and a : the potential increases 
when D increases whereas it decreases when a increases. The oscillatory behavior 
of the absorption cross-section, that follows as modes with higher angular mo- 
mentum numbers gradually come into dominance, is also affected by the presence 
of the Gauss-Bonnet term: the frequency of the oscillations increases with the 
Gauss-Bonnet coupling constant. 

3.2.2 Absorption cross-section for fermions 

For fermion fields, the ingoing flux is defined as the radial component of the 
conserved current = 2o~AB^f/2^i/2 evaluated over a two dimensional sphere 
(see f° r an explicit calculation). As in the case of scalar fields, this flux has 
to be calculated both at the black-hole event horizon and at spatial infinity in 
order to find \Aj\ 2 . Once this is done, the absorption probability for a mode with 
an angular momentum j is found to be 

\A 3 

Unlike the scalar case, we can not use the outgoing flux to determine the ab- 
sorption coefficient because we have only considered the upper component of the 
fermion fields. 

Figure H3 displays the numerical evaluation of the absorption cross-section for 
fermion particles on the brane. The high-energy limit of this cross-section is 



1 

B~i 



(19) 
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equal to the one for scalar particles and behaves in the same way with respect 
to the number of dimensions and to the Gauss-Bonnet coupling constant : the 
high-energy limit increases with the Gauss-Bonnet coupling constant whereas it 
decreases with the number of dimensions. When uth — > 0, the absorption cross- 
section tends to a non-vanishing value which behaves similarly to the Schwarzs- 
child case with respect to the number of dimensions : an increase in the value 
of D leads to an increase of the low-energy asymptotic value of the absorption 
cross-section (we refer the reader to U\ for an exhaustive discussion). The 
Gauss-Bonnet term leaves also a footprint on the asymptotic value of the cross- 
section when ujru — > : this asymptotic value decreases when the Gauss-Bonnet 
coupling increases. As in the high energy limit, the Gauss-Bonnet coupling con- 
stant changes the absorption cross-section in an opposite way than the number 
of dimensions does. 



3.2.3 Absorption cross-section for gauge bosons 

For gauge boson fields, the ingoing flux is derived from the trace of the energy- 
momentum tensor T^ w = 2o'ab (J A'B'^i B ^i B evaluated, as for fermions, over a 
two dimensional sphere. An explicit calculation of the flux at the black hole's 
event horizon and at spatial infinity leads to the absorption coefficient E] : 



' H 



1 



(20) 



Despite a few quantitative differences, we may clearly see from Fig. 01 that the 
absorption cross-section in the case of gauge bosons exhibits a similar behavior to 
fermions with respect to the number of dimensions and to the Gauss-Bonnet cou- 
pling constant. The high-energy regime asymptotic values are again suppressed 
with the number of dimensions and enhanced with the value of the Gauss-Bonnet 
coupling constant. The main difference occurs in the asymptotic regime of low 
energy: the absorption cross-section tends to zero when curjj — > 0, a behavior 
exhibited also in the pure Schwarzschild case [Oj. Although the asymptotic value 
of the absorption cross-section vanishes in this regime, the statement that an 
increase in the Gauss-Bonnet coupling constant leads to a decrease of the cross- 
section at low energy is still valid as it falls down to zero more rapidly for higher 
values of a. On the contrary, the cross-section approaches zero slower for higher 
values of the number of dimensions. 

From the previous results, we can summarize the effect of the Gauss-Bonnet 
term on the absorption cross-section as follows : for scalar particles, an increase of 
the Gauss-Bonnet coupling constant leads to an increase of the absorption cross- 
section over the whole energy spectrum. For particles with non-vanishing spin, 
an increase of the Gauss-Bonnet coupling constant leads to an increase of the 
absorption cross-section in the high-energy regime whereas it decreases it in the 
low-energy regime. Finally, the effect of the number of dimensions is the same as 
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Figure 4: Absorption cross-section for gauge bosons on the brane for Schwarzs- 
child- Gauss-Bonnet black holes as a function of the energy of the particle urn- 
Left Panel : the Gauss-Bonnet coupling constant a is fixed at 0.1 M~ 2 and the 
number of dimensions D at {5,6, 10}. Right Panel : the Gauss-Bonnet coupling 
constant takes the values {0, 0.1, 1, 10} M~ 2 and D is fixed at 6. 



for Schwarzschild black holes and affects the cross-section in exactly the opposite 
way that the Gauss-Bonnet coupling constant does. As a similar behavior is also 
exhibited by the temperature of the black hole in terms of its dependence on a 
and D, it becomes increasingly more interesting to compute the final emission 
rate to decide on the prevailing factor (temperature or absorption cross-section) 
that shapes the radiation spectrum for this particular black-hole background. 



3.3 Semi-analytical results for the high energy limit 

As it is known, and was numerically checked in the previous subsections, in 
the high-energy regime, the greybody factor assumes its geometrical optics limit 
value, which is independent of the spin of the emitted particle. This limiting value 
has been successfully used in both 4-dimensional and higher-dimensional black- 
hole backgrounds E] to describe the greybody factor of the corresponding 
black hole. This technique makes use of the geometry around the emitting (or 
absorbing) body and the same method can be applied for a Schwarzschild- Gauss- 
Bonnet projected black hole: for a massless particle in a circular orbit around a 
black hole, described by a line-element of the form ©, its equation of motion, 
P^Pn — 0, can be written in the form 

I ^L) 2 = 1 _ M^l (21) 

r dip/ b 2 r 2 ' 

where b is the ratio of the angular momentum of the particle over its linear momen- 
tum. The classically accessible regime is defined by the relation b < mm(r/\fh). 
The above result is valid for all spherically-symmetric Schwarzschild-like line- 
elements of the form Q with the structure of the particular gravitational back- 
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rr„ /( Trr'ir) 


a = 


a = 0.1 


a = 1 


a = 10 


D=5 


4.00 


4.59 


9.46 


51.2 


D=6 


3.07 


3.91 


7.64 


11.4 


D=7 


2.60 


3.53 


5.70 


6.62 


D=8 


2.31 


3.24 


4.53 


4.88 


D=9 


2.12 


3.00 


3.82 


3.98 


D=10 


1.98 


2.80 


3.34 


3.43 


D=ll 


1.87 


2.63 


3.00 


3.06 



Table 1: High-energy limits of the greybody factor on the brane in units of itr 2 H . 

ground entering the above relation through the metric function h(r). In order to 
compute the critical value of the parameter b, we must first minimize the function 
r/y/h] this takes place at r = r c , where r c is defined through the relation 



6Aira GM BH (D-3)(D-A) = 8tt(£> - l)GM BH 

(D-2)Q D ^ 2 r^ (D-2)fii,_ 2 rf-3' 1 J 

Then, b c is found by evaluating the function rj \/h(r) at the value r = r c . The 
corresponding area, a g = irb^, defines the absorptive area of the black hole at 
high energies and, thus, its greybody factor - being a constant, the emitting 
body is characterized at high energies by a blackbody radiation spectrum. Its 
exact expression comes out to be 



2 



rl L Q47raGM BH (D-3)(D-4:) 

2a(D -3)(D-A)\ l \ + (D - 2) Q D _ 2 rf" 1 



(23) 



For D = 4 and a = 0, r c = 3rn and the above expression reduces to the four- 
dimensional one, a g = 277rr^/4. For arbitrary values, however, of D and a, Eq. 
(1221) can not be solved in a closed form, and the value of r c can be found only 
numerically in terms of D, a, and, through Eq. (j3J), r#. In Table 1, we display the 
values of the greybody factor under the geometric optics approximation, in units 
of 7rr^, for various values of D and a, and for th = 1- It can easily been checked 
that these results are in excellent agreement with the high-energy behavior of the 
numerically computed greybody factors presented in the previous subsections. 



4 Radiation in the bulk 

4.1 Equation for scalar particles in the bulk 

Although all Standard Model (SM) particles are confined to our 4— dimensional 
brane, scalar particles which are neutral with respect to the SM charges are 
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allowed to propagate in the bulk. As a consequence, black holes can radiate spin- 
particles out of the brane and this radiation, although invisible to us, contributes 
to the total mass-loss rate. In order to make accurate estimates of the amount 
of energy radiated by the black hole on the brane, it is therefore necessary to 
determine the absorption cross-section and eventually the energy emission rate 
for those particles, too. 

The generalization of the Klein-Gordon equation to a D— dimensional curved 
space-time is given by \& ;/ ^ = 0, where /x runs over the D dimensions. In the 
spherically symmetric case, the following ansatz is used : 

(*, r, 9i, <p) = e-^P(r) Y e m (9 u <p) (24) 

where Y[ l {9 i ^ip) stands for the D— dimensional generalization of the spherical 
harmonics [32] • With this ansatz, the radial part of the Klein-Gordon equation, 
satisfied by the unknown function P(r), takes the form: 



W d (r°->h(r) d ~f] + 



h(r 

r D-2 c i r y dr / 



P = 0. (25) 



Apart from an implicit dependence upon the number of dimensions through the 
metric function h, the bulk equation exhibits also an explicit dependence through 
the contribution of the D — 2 compact dimensions entering the kinetic term and 
through the eigenvalues associated to the D— dimensional spherical harmonics. 
These contributions modify the exact form of the centrifugal potential, however, 
its general behaviour remains the same vanishing both at the black hole horizon 
and spatial infinity. Therefore the solution of the above equation at these two 
asymptotic regimes will again be given in terms of plane waves. 

Although the basis for the Hawking radiation process in the bulk is the same 
as for emission on the brane, both analytical and numerical studies of Eq. fl25(l 
have shown that varying the number of dimensions have different quantitative 
consequences on the bulk greybody factor and emission rate, for Schwarzschild 
[7] and Schwarzschild-de Sitter black holes |14j . As in those studies, we expect the 
number of dimensions as well as the Gauss-Bonnet coupling constant to strongly 
modify the value of the absorption cross-section for scalar particles in the case of 
a Schwarzschild- Gauss- Bonnet black hole, too. 



4.2 Numerical results for the absorption cross-section 

We numerically evaluate the absorption cross-section for emission in the bulk 
by following the same method as in the case of brane emission : Eq. ()25|) is 
numerically integrated starting from the vicinity of the black hole's event horizon, 
at r ini = r H +e, where the boundary conditions (|13j ) -(|14j ) have been again applied, 
and propagating the solution until spatial infinity. The analytical asymptotic 
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Normalized energy of the particle [a> r H ] 
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Figure 5: Absorption cross-section for scalar fields in the bulk for Schwarzschild- 
Gauss-Bonnet black holes as a function of the energy of the particle urn- Left 
Panel : the Gauss-Bonnet coupling constant a is fixed at 0.1 M~ 2 and the number 
of dimensions D at {5, 6, 10}. Right Panel: the Gauss-Bonnet coupling constant 
takes the values {0, 0.1, 1, 10} M~ 2 and the number of dimensions is fixed at 6. 



solutions at spatial infinity, given by j3] 

are then fitted on the numerical results and the amplitude of both ingoing and 
outgoing modes are extracted from the fit parameters. Once again, the absorption 
probability \Ae\ 2 is linked to these amplitudes via Eq. ()18j) . and the absorption 
cross-section is obtained using the D— dimensional generalization of the optical 
theorem given in Eq. (JHJ). 

Our results are shown in Fig. EJwhere the absorption cross-section is expressed 
in units of the black hole's area Ae>_2, i.e. the area of a (D — 2)— dimensional 
sphere with radius r# 

A D . 2 = r£" 2 (2tt) TT^ 3 )/ 2 r (^i) _1 . (27) 

These plots exhibit the same qualitative behavior as for the brane emission case : 
the absorption cross-section tends to a non-vanishing value when ur^ — ■> 0, in- 
creases with the Gauss-Bonnet coupling constant while it decreases for higher D 
values. From a quantitative point of view, the number of dimensions leaves a 
clear footprint on the asymptotic value in the low energy regime as a — > A0-2, 
which only depends on the size of the black hole and the number of dimensions, 
a behavior similar to the one found in the D— dimensional Schwarzschild black 
hole case [7j. The effect of the Gauss-Bonnet coupling constant, on the other 
hand, is greatly enhanced when compared to what happens in the case of brane 
emission. The significant enhancement of the bulk absorption cress-section, as 
the Gauss-Bonnet coupling constant increases, will play a significant role in the 
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(T An o 
u g 1 J ^U—Z 


a = 


a = 0.1 


a = 1 


a = 10 


D=5 


1.70 


2.09 


6.19 


77.7 


D=6 


1.77 


2.87 


11.0 


24.5 


D=7 


1.85 


3.97 


13.2 


19.1 


D=8 


1.93 


5.31 


14.6 


18.1 


D=9 


2.01 


6.81 


15.8 


18.3 


D=10 


2.08 


8.40 


17.0 


19.0 


D=ll 


2.16 


10.0 


18.3 


19.9 



Table 2: High-energy limit of the greybody factor in the bulk in units of Ad~2- 

determination of the bulk-to-brane emission rates ratio to be computed in Section 
5.3. Another novel characteristic, that arises in the case of a Schwarzschild-Gauss- 
Bonnet black hole, is the fact that the bulk absorption cross-section is enhanced 
by the number of extra dimensions from the intermediate energy-regime and 
onwards - in the case of a pure D— dimensional Schwarzschild black hole, the 
absorption cross-section in the bulk was heavily suppressed with the number of 
dimensions in the low- and intermediate-energy regimes allowing only for a small 
enhancement in the high-energy part of the spectrum [7| IS} I14j. 



4.3 Semi-analytical results for the high energy limit 

For emission in the bulk, the high-energy limit of the absorption cross-section 
is again given by the geometric optics approximation. The classically accessible 
regime for D— dimensional spacetimes of the form described in Eq. (0) is again 
given by b < mm(r/\/h) [HHj, with the minimum radial proximity b c obtained as 
in Section 3.3. What nevertheless changes, in the case of propagation in a higher- 
dimensional black-hole background, is the expression of the absorptive area of 
the black hole as a function of the critical radius b c : as it was pointed out in Ref. 
[7j, when carefully computed, this area comes out to be 

2n „(D— 4)/2 

p ~ (D-2) r[(D-2)/2] c ' 1 ' 

This area stands for the absorption cross-section a g , or greybody factor, for par- 
ticle emission in the bulk at high energies. 

In Table we give some indicative values of the greybody factor, in units of 
the area of the .D- dimensional horizon Ao-2-, for emission of scalar fields in the 
bulk at high-energies, for particular values of D and a, and for th — 1- Once 
again, these asymptotic values are in good agreement with the numerical results 
obtained in the ur^ — > oo limit, and reveal the enhancement of the absorption 
cross-section in this energy regime, in terms of both the Gauss-Bonnet coupling 
constant and the dimensionality of spacetime. 
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5 Particle emission spectrum 



In this section, the Hawking radiation spectrum emitted by a Schwarzschild- 
Gauss-Bonnet black hole is computed using the exact numerical results derived in 
the previous sections for the absorption cross-section. The spectrum is evaluated 
for scalar particles, emitted on the brane and in the bulk, and for fermion and 
gauge boson particles emitted on the brane. We also compare the number of 
scalar particles emitted on the brane and in the bulk as a function of the Gauss- 
Bonnet coupling constant, number of dimensions and mass of the evaporating 
black hole. It was pointed in Ref. [Hj that the calculation of the temperature of 
a Schwarzschild-de Sitter black hole is not simply given by the surface gravity as 
it is not possible to place an observer at spatial infinity. Such a modification in 
the temperature computation is a direct consequence of the intrinsic curvature of 
space-time. In the present case, the space-time is asymptotically flat and we can 
safely use the surface gravity to compute the Schwarzschild-Gauss-Bonnet black 
hole temperature given by equation P|). Consequently, the radiation spectra 
derived hereafter are relative to an observer infinitely far from the black hole. 

Up to this point, all thermodynamical quantities of the Schwarzschild-Gauss- 
Bonnet black hole, as well as the greybody factors, were conveniently expressed 
in terms of the event horizon radius r#. Comparing the values of these quantities, 
as the parameters D and a varied while r# remained fixed, was indeed an easy 
task. It also allowed for a direct comparison of our results to those derived 
previously in the literature under the same normalization. However, from an 
observational point of view, the horizon radius is not a convenient parameter to 
measure (and thus to fix), as different values of D and a correspond to black holes 
with different masses. Therefore, for observational convenience, in what follows, 
the various radiation spectra will be derived and compared at fixed black hole 
masses. 



5.1 Radiation spectrum on the brane 

For massless particles emitted on the brane, the flux spectrum, i.e. the number 
of particles emitted per unit time and unit frequency by the black hole, is given 
by Eq. (JJJ) evaluated for D = 4 : 

d 2 N s 1 4trx s H u 2 

dtdu (2vr) 3 e-/ T « - (-l) 2s ' 1 ' 

where s is the spin of the emitted field. Figures El El an d El depict flux spectra for 
scalars, fermions and gauge bosons, respectively, emitted by a higher-dimensional 
Schwarzschild-Gauss-Bonnet black hole. As the dependence of these spectra on 
the number of dimensions is similar to the one obtained in the case of a higher- 
dimensional Schwarzschild black hole, that has been exhaustively studied in the 
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Figure 6: Flux spectrum for scalar fields emitted on the brane by a Schwarzschild- 
Gauss-Bonnet black hole as a function of the energy u. The Gauss-Bonnet cou- 
pling constant takes the values {0, 0.1, 1, 10} M~ 2 and the number of dimensions 
is fixed at 6. The mass of the black hole is taken equal to Mbh = 10 M* 
(upper left), M BH = 100 M* (upper right), M BH = 1000 M* (lower left) and 
M BH = 10000 M* (lower right). 



literature already El E|, here we choose to focus on the effect of the Gauss- 
Bonnet coupling constant instead. We thus fix the number of dimensions at 
D — 6, and assign to a the values {0, 0.1, 1, 10} M~ 2 . In the case of scalar fields, 
shown in Fig. 6, the particle emission spectrum was also derived for four different 
values of the black hole mass, i.e. for {10, 100, 1000, 10000} M*. The spectra for 
fermions and gauge bosons are shown only for the two extreme mass values as 
their behavior is qualitatively identical to the one seen in the scalar case. 

From these plots, we may clearly see the non-monotonic effect that the Gauss- 
Bonnet coupling constant has on the number of emitted particles. We can also 
clearly see a strong dependence of the spectrum on the mass of the evaporating 
black hole. A number of competing factors contribute to the non-monotonic 
behaviour in terms of a displayed in Fig. El although the greybody factor for 
scalar fields is greatly enhanced in units of the horizon area, the horizon radius 
itself shrinks under an increase in the value of the Gauss-Bonnet coupling, when 
we keep the mass of the black hole fixed; in addition, the temperature of the 
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Figure 7: Flux spectrum for fermion fields emitted on the brane by a Schwarzs- 
child- Gauss-Bonnet black hole as a function of the energy uj. The Gauss-Bonnet 
coupling constant takes the values {0, 0.1, 1, 10} M~ 2 and the number of dimen- 
sions is fixed at 6. The mass of the black hole is Mbh = 10 M* on the left panel, 
and M BH = 10000 M* on the right panel. 



black hole as a function of a is also non-monotonic - a simple numerical analysis 
shows that it initially decreases with a but, after a turning point is reached, it 
starts increasing again. All the above factors tend to compensate each other, 
with the outcome of this process greatly depending on the black hole mass value: 
for light black holes, the reduction of the horizon radius, as a increases, is a 
significant one which, in conjunction with the decreasing temperature, leads to 
the suppression of the number of particles emitted - only after the turning point 
of the temperature is reached does the spectrum show some enhancement. For 
heavier black holes, the effect on the value of the horizon radius, as well as the 
decline of the temperature with a, tends to become less and less important; the 
increase of the absorption cross-section then prevails over all other factors, thus 
leading to an enhancement of the radiation spectrum as the value of a increases. 

Similar arguments hold for the case of the emission on the brane of fermions 
and gauge bosons. In both cases, for low-mass black holes, the decrease of the 
horizon radius and of the temperature leads to a suppression of the number of 
particles emitted, even in the intermediate and high-energy regime where the 
greybody factor was greatly enhanced as a function of a; only for large values of 
the Gauss-Bonnet coupling constant does the spectrum shows signs of enhance- 
ment as a reflection again of the change in the temperature behavior. For black 
holes with larger masses, the enhancement of the greybody factor with a, found 
for the intermediate and high-energy regime, prevails again and eventually leads 
to a small enhancement in the number of particles emitted. This enhancement 
is slightly smaller than the one found in the case of the scalar fields as a result 
of the radically different behaviour of the greybody factor for different fields at 
the low-energy regime: while the enhancement, in terms of a, seen in this energy 
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Figure 8: Flux spectrum of gauge bosons emitted on the brane by a Schwarzschild- 
Gauss-Bonnet black hole as a function of the energy u. The Gauss-Bonnet cou- 
pling constant takes the values {0, 0.1, 1, 10} M~ 2 and the number of dimensions 
is fixed at 6. The mass of the black hole is Mbh = 10 M* on the left panel, and 
M B h = 10000 M* on the right panel. 



regime gives a boost in the spectrum of scalar fields, its suppression for fermion 
and gauge bosons restricts the enhancement of the flux spectrum at higher energy 
regimes. 

Let us finally note here that the non-monotonic behavior of the temperature 
of a Schwarzschild-Gauss-Bonnet black hole as a function of the Gauss-Bonnet 
coupling constant is clearly reflected on the value of the energy peaks in the 
particle spectra. For instance, the temperature of a 6— dimensional black hole 
with a mass of 10 M*, is roughly equal to Tjj = {0.14, 0.11, 0.06, 0.5} M* for 
a = {0,0.1, 1, 10} M~ 2 . Focusing, for example, on the scalar particles spectrum 
emitted by a black hole with a mass of 10 M* (upper left panel of Fig. EJ), we see 
that the peaks of the emission curves shift towards lower or higher energy values 
accordingly. Moreover, they correspond to energies lying in the range ~ 1.4 x T# 
to ~ 2 x Th-i depending on the exact value of a, which is in good agreement with 
the prediction for the location of the energy peak for a pure black-body spectrum 
at ~ 1.6 x T H . 



5.2 Radiation spectrum in the bulk 

For massless particles emitted in the bulk, the Hawking radiation spectrum is 
given by : 

d 2 N 1 Qfl-2 a(uj) u) D ~ 2 

dtdu (2tt) d - 1 e"' T «-l ' 1 ' 

where the spin index is omitted as only scalar fields can propagate into the bulk 
(gravitons are not considered in this work). By using our numerical results for 
the bulk absorption cross-section, derived in section 4.2, we may now proceed to 
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Energy of the particle a> in Planck units Energy of the particle co in Planck units 

Figure 9: Flux spectrum for scalar fields emitted in the bulk by a Schwarzschild- 
Gauss-Bonnet black hole as a function of the energy u. The Gauss-Bonnet cou- 
pling constant takes the values {0, 0.1, 1, 10} M~ 2 and the number of dimensions 
is fixed at 6. The mass of the black hole is Mbh = 10 M* on the left panel, and 
M B h = 10000 M* on the right panel. 



derive the particle emission spectrum for emission in the bulk of scalar fields. The 
absorption cross-section will be expressed in units of the D— dimensional black 
hole area Ad~2 = but, as in the previous subsection, we choose the 

mass of the black hole as its characteristic parameter that remains fixed. 

The flux spectra for scalar particles emitted in the bulk by a 6— dimensional 
black hole, with a mass of 10 M* and 10000 M*, are shown in the two pan- 
els of Fig. for a set of values of the Gauss-Bonnet coupling constant given by 
{0, 0.1, 1, 10} M~ 2 . The spectra exhibit the same characteristics, in terms of their 
behavior as functions of the Gauss-Bonnet coupling constant, as in the case of 
scalar fields emitted on the brane. When the mass of the black hole is only a few 
times above the Planck scale, the decrease in the horizon value and the temper- 
ature of the black hole, when a increases, eventually leads to the suppression of 
the flux spectrum, despite the enhancement of the absorption cross-section over 
the whole energy regime in units of the black hole's area. For significantly heavier 
black holes though, the magnitude of this enhancement overcomes the decrease 
in the area and temperature of the black hole, leading to a significant increase in 
the number of particles emitted, especially for high values of a. 

5.3 Bulk-to-brane ratio 

The bulk-to-brane relative emission is of great importance for experimental inves- 
tigations: particles emitted in the bulk cannot be detected, and the correspond- 
ing energy will be permanently lost. It was shown in Refs. [HE] that this ratio 
remains lower than one for D— dimensional Schwarzschild and D— dimensional 
Schwarzschild-de Sitter black holes. Such results clearly favour the scenario of the 
experimental detection of small black holes through the emission of Hawking ra- 



21 





Figure 10: Bulk-to-brane relative emission rate as a function of the energy to. 
The Gauss-Bonnet coupling constant takes the values {0, 0.1, 1, 10} M~ 2 and the 
number of dimensions is fixed at 6. The mass of the black hole is taken to be 
M B h = 10 M* (upper left), M BH = 100 M* (upper right), M BH = 1000 M* (lower 
left) and M BH = 10000 M* (lower right). 



diation since most of their energy will be spent on the emission of brane-localised 
particles. However, as it was mentioned in Section 4.2, for Schwarzschild-Gauss- 
Bonnet black holes, the effect of the Gauss-Bonnet coupling constant on the 
absorption cross-section is more important for particles emitted in the bulk than 
for particles emitted on the brane. Therefore, we can not exclude the possibility 
that such black holes can decay mainly in the bulk channels instead of the brane 
channels. We expect this to happen especially in the high-energy regime, where 
the cross-section is more significantly enhanced, and for heavy black holes, for 
which the effect of the greybody factor dominates the spectrum. 

In Fig. El we present the outcome of our attempt to check this possibility 
numerically by using our results for the bulk and brane particle emission rates 
for scalar fields. The four panels of this figure show the bulk-to-brane ratio for a 
6— dimensional black hole with a mass equal to Mbh — {10, 100, 1000, 10000} M*, 
respectively, and for different values of the Gauss-Bonnet coupling constant. As 
we see, in almost all cases, the relative emission rate remains lower than one, 
although it asymptotically approaches this value in the high-energy regime. As 
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Figure 11: Bulk-to-brane relative emission rate as a function of the unidimen- 
sional parameter ujth for 6— dimensional black holes with a given size. Left Panel : 
the Gauss-Bonnet coupling constant is taken at {0, 0.1, 1, 10} M~ 2 and the num- 
ber of dimensions is fixed at D — 6. Right Panel : the Gauss-Bonnet coupling 
constant is fixed at 0.1 M~ 2 and the number of dimensions is taken at {5, 6, 10}. 



expected, this ratio exceeds unity only in the case of a black hole with mass of 
10000 M* and only in the high-energy part of the spectrum. Therefore, it seems 
highly possible that, as we increase further the value of the black hole mass, and 
for sufficiently high Gauss-Bonnet coupling constant, the bulk channel becomes 
increasingly more dominant, leading eventually to a significant loss of the black 
hole energy into the bulk. 

An upper limit on the bulk-to-brane relative emission rate can be obtained 
by looking at the spectrum of black holes with fixed horizon radius instead of 
a fixed mass. Figure ITTT left panel) shows the spectrum of such a black hole. 
Since r# is now kept fixed, the different curves shown in this figure for different 
values of the Gauss-Bonnet coupling constant correspond to black holes with 
different masses, according to Eq. Q. However, there is one case where all 
these curves would describe a black hole with the same mass, and that is the 
case of a black hole with an infinite mass - in that case, any change in a would 
produce a negligible change in the black hole mass. Therefore, Fig. ^2 can be 
used to derive an upper limit for the bulk-to-brane ratio. As we may see, for 
high values of the Gauss-Bonnet coupling constant, the relative emission rate 
becomes indeed greater than one. This value is attained over the whole high- 
energy regime, and can therefore lead to a major loss of the black hole energy 
into the bulk (or, even the decay of the black hole mainly through the bulk 
channel) via the emission of particles with uj > l/r#. We should, however, note 
here that such a scenario is not experimentally favored as it requires very high 
values for the Gauss-Bonnet coupling constant and black hole masses above the 
values that could be reached by the next-generation colliders. The behavior of 
the bulk-to-brane relative emission rate for an intermediate value of the Gauss- 
Bonnet coupling constant, i.e. a = 0.1 M~ 2 , is shown in Fig. ITTT right panel) as 
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a function of the number of dimensions : as in the case of a pure D— dimensional 
Schwarzschild black hole [7j it remains lower than one and decreases for higher 
dimensional black holes. 

6 Conclusions 

Looking for gravitational effects beyond general relativity is a major experimental 
challenge. Lovelock gravity (especially its curvature-squared Gauss-Bonnet term) 
is an interesting way of extending the Einstein-Hilbert action in a string-inspired 
approach. As was shown in [T^], Schwarzschild-Gauss-Bonnet black holes, that 
could also be formed at colliders, would allow an accurate investigation of the 
associated coupling constant. This paper goes beyond the approximation made 
in ^1 for the expression of the Hawking radiation emission rate of these black 
holes, and presents the complete, exact particle spectrum by including the missing 
material, i.e. the greybody factors. 

More analytically, in this paper, we have investigated the decay of a D— dimen- 
sional Schwarzschild-Gauss-Bonnet black hole via the emission of Hawking radia- 
tion. For this purpose, we have computed the exact expressions for the greybody 
factors for scalar, fermion and gauge boson fields emitted on the brane as well as 
for scalar fields emitted in the bulk by using numerical analysis. These results 
have revealed that the absorption cross-sections depend strongly not only on the 
energy and spin of the emitted particle and the dimensionality of spacetime but 
also on the value of the Gauss-Bonnet coupling constant. Interestingly enough, it 
turns out that the effect of this coupling constant on the absorption cross-section 
is exactly the opposite of the one of the dimensionality of spacetime : this quan- 
tity is enhanced over the whole energy regime for scalar fields (emitted both on 
the brane and in the bulk) while, for fermions and gauge fields, it is enhanced 
in the intermediate- and high-energy regime and suppressed in the low-energy 
regime. 

With these results at hand, we proceeded to evaluate the particle emission 
spectra for all types of particles. At this point, we switched from the convenient 
assumption of a black hole with a fixed horizon radius to the phenomenologi- 
cally realistic one of a black hole with a fixed mass - that allowed us to reveal 
and emphasize the importance of the mass of the black hole on the behavior of 
the radiation spectra with respect to the Gauss-Bonnet coupling constant. Our 
results have shown that, for light black holes, an increase in the value of the 
Gauss-Bonnet coupling constant leads to the suppression of the number of parti- 
cles emitted by the black hole per unit time and unit frequency, for all types of 
particles. This is the result of the ensuing decrease of the black hole's area and 
temperature dominating over any enhancement exhibited by the greybody factor. 
However, for heavy black holes, it is the enhancement of the greybody factor that 
overcomes the small, in this case, decrease of the black hole's area and temper- 
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ature, thus leading to a spectrum that is also predominantly enhanced with the 
Gauss-Bonnet coupling constant. 

The bulk-to-brane relative emission rates have also been computed for the 
emission of scalar fields that are allowed to propagate both on and off the brane. 
They were shown to be also black-hole-mass and coupling-constant-dependent. 
Whereas, for small or zero Gauss-Bonnet coupling constant, this ratio remains al- 
ways smaller than unity - thus favoring the detection of mini black holes through 
their decay mainly on the brane where our detectors are located - this situation 
may easily change for heavy black holes and sufficiently high Gauss-Bonnet cou- 
pling constant. Nevertheless, any significant effect of this type would demand 
black hole masses beyond the order of 10000 M* and a Gauss-Bonnet coupling 
constant of order 10 M~ 2 . Therefore, for black holes with masses attainable at 
next-generation colliders, their decay will take place mainly through brane chan- 
nels, even if quadratic curvature terms are taken into account in the theory. 

To conclude, let us note that the above results have been derived under the as- 
sumption that the extra spacelike dimensions can be considered as non-compact. 
This assumption might be violated if the black hole radius becomes comparable 
to or larger than the size of the extra dimensions. Depending on the number of 
dimensions, the latter varies in the range of a few fermis to a few millimeters 
in order to obtain a Planck scale of order a TeV without being in contradiction 
with any experimental observations. As black holes produced at next-generation 
colliders will not be heavier than 10 M*, their sizes will remain much lower than 
a fermi, and our calculations can therefore be safely used to compute the de- 
cay of higher-dimensional Schwarzschild-Gauss-Bonnet black holes. Finally, our 
results apply for particles with wavelengths lower than the size of the extra di- 
mensions, which then defines a lower limit on the particle's energy. Therefore, 
the interpretation of our results found for the very low-energy part of the derived 
spectra must be interpreted with care. This constraint however does not dimin- 
ish the importance of our results or the possibility of detecting the presence of 
quadratic gravitational terms in the theory: the effect of the Gauss-Bonnet term 
becomes manifest in a clear way over the whole radiation spectrum, especially in 
the intermediate and high-energy regime. 
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